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ABSTRACT 

A shear flow aerodynamic theory for steady incompressible flows 
is presented for both the lifting and non^lifting problems. The unique 
feature of the present theory is the consideration of the slow variation 
of the boundary layer thickness. The slowly varying behavior is treated 
by using the method of multi-time sciAes. The analysis begins with 
the elementary wavy wall problem rnd, through Fourier superpositions 
over the wave number space » the shear flow equivalents to the aero- 
dynamic transfer functions of classical potential flow are obtained. 

The aerodynamic transfer functions provide integral equations which 
relate the wall pressure and the upwash. Computational results are 
presented for the pressure distribution, the lift coefficient, and 
the center of pressure travel along a two dimensional flat plate in a 
shear flow. The aerodynamic load is decreased by the shear layer, 
compared to the potential flow, while the variable thickness shear 
layer decreases it less than the uniform thickness shear layer based 
upon equal maximum shear layer thicknesses. 
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perturbed characteristic length in aean flow direction 
Aerodynamic influence functions in physical space 
Aerodynamic influence functions in wave number space 

ft 

See Equation (16) 

Wing span 
Wing chord 

Surface Contour (Figure 1) 

Wavy wall amplitude 

Bessel fimction of second kind of order v 
Aerodynamic kernel functions in physical space 
Aerodynamic kernel functions in wave number space 
Characteristic length in mean flow direction 
See Equations (16) and (21) 

Exponent in shear layer velocity profile (see Eq.(4 )). 
Fluid perturbation pressure 
Wavy wall fluid pressure 
Wall pressure 

Wavy wall fluid pressure at wall 
H (o2 + y2)^^^ 
s VH 

Fluid perturbation velocities in x» y and z directions 
Shear flow mean velocity 
Free stream velocity 
Upwash s 

Wavy wall upwash 
See Equation (40) 
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X, y, * 

- Coorainate axes 

X', y'. *• 

- (see Figure 1) 


- See Equation (8) 
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- Fluid density 
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. Gaaaa function 
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. Botmdaxy layer thickness 

^nax 

- Maxiitum boundary layer thickness along chord 
(trailing edge) 

«» Y 

- Wave numbers in x and y directions 
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*«ax/^ 

V 

-i 1,1 in Bain text 


> fluid kinematic viscosity in appendix. 

Cl 

- 5 total lift / Cy 0 U^.c) 

C(. 

- 2 moment about leading edge / t y P • c^) 



I. INTRODUCTION 


To account for the presence of the boundary layer adjacent to any 
solid surface in a fluid flow theoretical analyses of the shear flow 
effect have been made by many researchers [1-20]. Conmon to all the 
above analyses is the assumption of the uniform thickness of the 
boundary layer initially present before the solid boundary surface is 
deformed. However, the boundary layer normally grows in the mean flow 
direction and one naturally would like to know how the boundary layer 
thickness variation affects the aerodynamic load on the solid surface. 

This paper deals with this problem but is restricted to the steady, 
incoapi*essible flow with a slowly varying boundary layer. Similar reasoning 
can be extended to unsteady and compressible flows. 

Since this report is an extension of Ventres' work [19], we first 
briefly review and quote his uniform thickness results and then develop 
a theory for the slowly varying boundary layer problem. 
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II. BRIEF REVIEW OF VENTRES' RESULTS 

Consider s steady* incompressible shear flow over a surface whose 
deflection is given as x ■ f(x,y). See Figure 1. The surface deflection 
produces a small disturbance from an initially parallell shear flow 
u*U(z)*v«w«0. The fimction U(x) is constant for z > 5 so that 
the shear layer is limited to the region 0 < z < 5 adjacent to the 
surface. 

With the choice of the mean flow U(z) Ventres was 

able to relate semr-analytically the upwash W(x,y) everywhere in the 
z ■ 0 plane to the perturbation pressure distribution on the wing p(x,y) 
(recall that we are dealing with the lifting problem) according to 


W(x>y) “ // 

Uj wing 


y-n) 


where W(x,y) s 3f(x,y) and fCx.y) is the wall deflection. The 


two-dimensional kernel function ^(x)* which is of fundamental importance, 
is siiown in Figure 2 for N • 7 and 11. Also shown is the potential flow 
result* 

KflcX- i 

ifX 

which is labelled in the Figure as N - All three curves have a common 
asymptote as x/^ ■* » * and in fact are essentially identical for x/^>2 
or so. Since K physically is the upwash caused by an Impulse pressure, 
this implies that the influence region of the shear layer is limited to 
a distance on each side of the source point comparable to the shear layer 






thickness, or the shear layer effect is "nearly local". This points 
up the possibility of accounting approximately for the effect of 
slowly varying boundary layer thickness by inserting a variable <(x) 
directly into Ventres uniform thicknesf shear layer theory. The 
problem is tiien how to insert the variable 6(xl appropriately. The 
rest of this paper will demonstrate the reasonable way of doing it 
using the multi-time scale concept [21]. 
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HI. SHEAR FLOW WITH SLOWLY VARYING BOUNDARY UYER THICKNESS 


(ft) Governing Equations And Boundary Conditions 

Detailed order of nagnitude analysis given in the Appendix shows 
that the perturbation pressure p(x,y«z) satifies 

V2p - 2 dz ^ - 0 (1) 

U 3z 

The associated boundary conditions are [12, 19} 


3z 


• -oU (z = z.) 3w 
3x 


( 2 ) 


on the solid surface as z^ 0, where 

I . 

Iz^Zq 3X 
and the finiteness condition 

p -► 0 as z • (3) 

We assume the mean flow. 



U 




1 L ^(x.y) 

2 1 5(x,y) 


(4) 


While assuming a discontinuity in the mean flow velocit)^ gradient across 
z ■ 4(x,y), we impose the continuity of the pressure and the pressure gradient 


across the boundary layer edge* i.e 


at I - 6(x*y) 

pCx,y,z » «‘Cx.y)) - p(x,y,z ««*(x,y)) 

«nd^£ (x,y.z »«"(Xiy)) ■ a£Cx,y,z»fi*Cx,y)) (S) 

3z 3z 

Froa Equations(l) and (4) the pressure within the shear layer 
region satisfies 


^ ^ . 0 ( 6 ) 

Nz 3z 

which will be solved along with the wall boundary condition. 

Equation (2), i.e. 


i£ 

3z 


- P U2 ill. 


as Zq 0 


( 7 ) 


and another boundary condition at the boundary layer edge to be derived 
from the field equation outside the boundary layer satisfying the 
finiteness condition far away from the solid surface (see Equation (14)}. 
(b) Normalization 

In view of the slowly varying behavior it is convenient to write 
the boundary layer thickness 
« - 6 (eC, e’n) 

idiere e << 1 and e* « 1. Here, e and c' are non-dimensional para 
meters characterizing the slow variation of the boundary layer thickness 


of the 

OKlGiiv.J. ^ 


- ^ - 


in the X nnd y directions respectively, e,g. e « and c* * 

i i» th« trailing edge shear layer thickness which is 

aax/^ max 

usually the maxistus along the chord of an airfoil and b is the wing 
span. C and n are non-dinensional coordinates defined as 


C = X , n • 31 

y s 

■ax max 


C 


TCeC»e*n) 


( 8 ) 


The definition for t has the advantage that all toundary conditions 
are applied on constant values of c , e.g. c 0, ■ 1 and •• 
(corresponding to z « 0, 6(x,y), and 
(c) Wavy Wall Problem • Elementary Solution 

Consider a wavy wall whose profile is describ<d by the real 
part of the complex function 


f Cx.y) . ? e 


? 9 i(«€+yn) 


(9) 


where o ■ o d. 


. y " y K 


max ' ' max. 

The wavy wall will generate a perturbation pressure field of the following 


form. 


P * p(C,e5.c'n; o,y;6»e') e 


iCoc+vn) 


( 10 ) 


in which the slight amplitude modulation due to the slow boundary 
layer variation is manifested by the functional dependence of p on 
e{,c*n e ^ e'. We further assume the following series expansion 

f or p since c « 1 and c* « 1[21]. 


I 


♦ epjCc,c£»e*n; o,y^ 


(U) 


♦ 6'Pj*(;»cC,£Vi;a,^ ♦ H.O.T. 

Rewriting Equation (6) in terns of C* C rather than x, y» z by 
means of Equation (8) and substituting the assumed wavy wall pressure 
solution. Equations (10) (11). into the resultant equation, one can 

show that the lowest order governing equation for is 


»Fo . 2 . fs \ 

^ Xt 35 J 


(12) 


where 4 is •‘h*' variable boundary layer thickness. The higher order 
equations, which are not shown here, have lower order solutions as 
their forcing functions. 

The required wall boundary condition for solving Equation (12) 
is obtained by substituting Equations (9) and (10) and (11) into 
Eo;iation (7) and identifying the zeroth order relation, namely 

-2 _ V / 4 


3-5 


p/o 0 v UV 0 

( J ) vr- 7 , 

z»Zq max ' max 




(13) 


as Cq 0. The higher order wall boundary conditions are trivial. 




0 for all i i 1 . " ^0 ° 


C-C, 


Another boundary condition at the edge of the shear layer for solving 
Equation (12) is obtained by setting N • «> in Equation (12) and applying 


tha finiteness condition a-t c 


• S 


i.g. 



( 14 ) 


Everywhere wisidt the boundary layer and in particular at the boundary layer 
edge ? = I because of the continuity of the pressure and the pressure 

gradient across c • 1. 

Equation (12) subject to boundary conditions (13) and (14) gives 
the wall pressure solution [19] » 


5«. ca.i.su.y)} • AC<.,Y.«U.y» 

“i \^7 


(15) 


where 


A ■ i 


a 

R 


2/jj 

^ J ^ . L(Rfi) 
R6 


A^- vr(i“v) 


r(i+v) 




, c® tW 


L(R«) • 


(16) 


W*(g,Y) » i o *(2ir)^ 

‘‘^1 

The solution is the same as Ventres’ except his uniform thickness has 
been replaced by the variable thickness 5(x,y). Henceforward, the 0 
subscript on p is dropped, for simplicity. 


- » - 

Thus we h«ve obtained the pressure aaplitude function for the 
wavy wall problim using the shear flow Dodel of slowly varying boundary 
layer. The pressure aaq^litude is proportional to the upwash ^ 
as expected in a linear theory. 

For an arbitrary wall deflection we can use the superposition 
property of a linear theory to form Fourier integrals. However the 
Fourier superposition should be performed carefully as shown in the 
following two sections. 

(d) Pressure Load for an Arbitrary Wall Deflection - Non>lifting Problem 

We have just derived the lilnoar relation for wall pressure 
amplitude (IS) for a wavy wall problem. Through the linear super- 
position of the wavy wall solutions over all possible wavr^. numbers, 
the wall pressure due to any properly behaved wall deflection is 

m m 

Cx.y,«Cx,y)) - / / p^ (o,Y;6Cx,y)) dody (17) 

5Uj2 PU^ 

Substitution of Equation ( 1 £) into Equation ( 17 ) yields 

P, tx.y,4(x,y))=' J.y ICo,y:S(x,yJ) If Csi.y) . £18) 

“l 


To write Equation (18) as on integral in physical space, we define 
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AU'.y'*«0t.y))s 


which iiqplies 



^-ifr.x'+Yy') 

0 


aady 


(18a) 


A(a,y;fi(x,y)) - // ACx',y';« Cx»y)) ^ dx’dy' 

In additioa, we have 

W (x.y) « V // i![lC“»Y) dadY 

“i yy - 

Using (18a) and (18c), one can easily show that Equation (18) is 
equivalent to 


(18b) 


(18c) 


*^w (x,y,5(x,y)) = // A(x-x*, y-y’;8(x,y)) W (x'»y*) dx*dy' (19) 

- Uj 

Altoxnatively, one can arrive at Equation (19) by using the Fourier 

convolution theorem directly from Equation (18) with the recognition 

that Equation (18a) is a definition af A(x',y*j <S(x,y)) in terms of 

A(e,Y»8(x,y)) and 8(x,y) is treated as irrelevant to the Fourier integrals. 

Either of Equation (l8) and Equation (19) can be used to calculate 

the wall perturbation pressure for a given upwash W(x,y). However, a 

drastic reduction of computation time can be made by recognizing that 

the quantity (1 - £(x,y))-is everywhere smaller than one except near 
5 

max 

the leading edge where 8(x,y) -*■ 0. Hence, we can <,'xpand the function A 

in Equation (18) in terms of a power series of (1 - 6_(x,y)); the 

''max 

resultant pressure solution presumably will be valid everywhere except 
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REPROBUC 


- in’rr / .t-. ttytt. 

-t wr THMJ 
■ ^ io Poor 


mar the leading edge where our theory is not expected to be accurate 
in any event because of the rapidly changing boundary layer thickness 
It can be shown that 




( 20 ) 


where 


^0* 


A * 
^1 


(*w) 


2 / 


N 


>• »w 


2 / 

*v»te) b 1 .axj 


■•»W ■ ‘v«W ♦ * v-l»W 




-y ^ 1-v^^max^ 

).. ^COS./^ /- 


♦ I ,.<M, 




*14 

) 

m 


(21) 


Substitution of Equation (20) into Equation (18) and the use of the 
Fourier convolution theorem gives the wall pressure. 


(x*y) =■ P„o ♦ A-« C^X P^I (x,y)*0A-l(x^’j^ 

t0fj2 ' ^max' pup '^max ' 


( 22 ) 


where 



Aq and Aj are the Fourier inversion of and Aj^* defined in Equation 
(21) and can be evaluated nuaerically, i.e. 

Aj(x,y) . ( I 1/ V * i<«'Y 

AiCx.y) • (l // A,* . d=dv 

Here, A^ is the sane as Ventres' result except bis uniform thickness 

has been replaced by the maximum thickness 5 A, is a new function 

used to calculate the pressure correction due to the slowly variation 

of the boundary layer, but it is independent of the detailed variation 

of the boundary layer which enters the pressure solution in the factor 

(1- in Equation (22). One can show that A. vanishes as 6 

®nax ^ 

approaches zero, a natural result because of the loss of the "variation" 
of the boundary layer as the boundary layer vanishes. 

Thus we have derived a schons for calculating the pressure load 
on any solid surface which is deformed slightly from its flat position 
and exposed to a shear flow whose boundary layer thickness slowly 
varies. In aeronautical jargon, we have solved the non-lifting problem 
because the equations we have found (Hqxiations (19) and (23)) require 
the upwash W(x,y) being given everywhere in the z =* o plane. For the 
technologically more important lifting problem the above Fourier super- 
position of the pressure load alone is not sufficient and a slightly 
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different procedure is needed as will be seen in the next section;. 

(e) The Pressure Load for an Arbitrar/ Wall Deflection > Lifting Problem 

For the lifting problffis, it is convenient to write Equation (IS) as 

( i f Pw C«,rHOt.y» C24) 

“i 

where Ic is the reciprocal of X defined in Equation (16) and is the same as 
Ventres' kernel with his uniform thickness replaced by the variable 
«(x,y), i.X^' 

K - Ay R /R« N^^y'LCRd) (2S) 

r a ^ 2 

Multiplied by o^C“^*Yy) integrated over (o.y) space. Equation (24) 
results in the following upwash equation. 


w Cx>y) «/i \^ // e ^-Ow+vy) ^ 26 ) 

^ Uj \?ff/ Uj 

- // k (o ,Yi 4 (x,yXX ?„(£»♦ Y;<S(x,y))l dady 

k pup 

In addition, the superposition result Equation Cl7) still holds for 

the lifting problem, i. e. 

I 
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Which can be formally inverted to give 


p^(a.y;«Cx,y)) 



Cx*,y';«Cx.y)) 

2 


^-iCax’ny’) 


dx’dy' 


(28) 


In obtaining Equation (2^)* fiCx*y) has been treated as irrelevant to 
the integral operation. Note that p^^Cx* ,y' Cxiy)) is not the real 
wall pressure as it should have been if x'« y* are replaced by x« y. 
Substituion of Equation (28X into Equation (26) gives the familiar kernel 
fiutetion form for the lifting problem. 


n (X,y) « // K(x-x'. y-y';«Cx,y)) p„ Cx’,y' ;6(x,y)) dx'dy' 


«i 

where 


Wing 




K(x,y,6(x,y))= (l\^ !! KCo,Y;6Cx,y))e"^^“*'*‘^^^dadY 
V 2Tf} -« 


(29) 


(30) 


The domain of integration in Equation (29) is within the wing surface 
area because p^ = 0 everyv*here off the wing for the lifting problem. 

The kernel function inversion. Equation (30), can be evaluated semi-ana- 
lytically according to Ventres . A convenient splitting of K 
into two parts is 


Ki>K, 


C3D 


where 


v2/ 


A^R/R6y'N(^ -n 
V LCR6) ' 



can be inverted numericaXiy to give and K 2 can be inverted 
analytically to give so that 

K - Kj ♦ K 2 (32) 

For two dinensional shear flow, 

K 2 Cx, 6 (x)) = a / 4 Y^N rC2v)eosjir 

and for three dimensional flow 

2/ X 

K,(x.y,S(x.y)) = - 'n / U itl 

Siailar to the non-lifting problem an expansion of the tranformed 

kernel Ic in terms of the quantity (1 - 6 (x,y) can be made. 

5 

max 





where 

r" a 2 ' 


(33) 


( 34 ) 
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V (» <«ax) 


l''" /Wr,ax> 

Vl- )/ 


Substitution of Equation (33) into Equation (30) gives 


K(x,y, (x,y)) » K^(x,y) / i-g Cx,y)^ K^(x,y)»on-6 Q^yl y 


aax 


max 


tAiare and Kg are the Fourier inversion of K^* and Kg* 
respectively, i.e. 

s &>•)') ■ /iV 

l 2irj -m 
and 


KgCx.y) =(l V // Kg* e ^t«x*yy) 


C3S) 


(36) 


K^ is the same as Ventres' K with his uniform thickness replaced by 

5 . The numerical inversion of K_* to give K« is straight forward, 

max e e 

Both K, and K„ for two dimensional shear flow are shown in Figure 3. 
A o 

We fiurther assume that 


P^, (x,y;3(x,y) = p^q (x.y) + / I 


iCx.yK (x,y)»0a ~3(x. y)\^ (37) 
^max ^ raax ^ 


which is presumably allowed as far as S(x,y) is not anywhere near zero. 

Substituting Equations (33) and (371 into Equation (29)., one can equate 

the terms of zeroth and first powers of /U - ^ ^*^?\ to yield the following 

' °max ' 


two kernel function forms; 


wing 

surface 


_i!o 


2 


and 

w flCx.y) - // y-y’) Pw Cx'.y') dx*dy» (39) 

_ wing ''l 

Uj* stirface pU^ 


where is an equivalent upwash in tenas of the lowest order solution 
P„ (x,y) i.o. 


W (x,y) = // - K.(x-x», y-y*) (x',y*) dx'dy' (40) 

Wing area “ 

“l ^ 

The dojaains of integration of the integrals in Equations (39) and 

(40) cover only the wing surface area because the pressure is zero 
everywhere off the wing surface for the lifting surface problem. 

Then the solution procedures are 

(1) to solve for p, (x,y) from Equation (38) for a given upwash 
W(x.y), 

(2) to calculate the equivalent upwash W (x»y) from Equation (40), 

(3) to solve for p (x,y) from Equation (59) , and 

'^1 

(4) to fozm the final solution using Equation (37) . 

Thus a lifting surface theory for variable thickness shear flow 
is con^jleted. It is seen from Equation (37) that the local pressure 
depends solely on its local boundary layer thickness for a slowly 
varying boundar)' flow to the lowest order approximation. However, 
Equation (29) says that not only the variable 5(x,y) should be 
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included in the kernel function K but also the wall pressure 
appears in the integral equation as the fictitious ^essure 
5(x,y) which is not the trie wall pressure P^^(x,y; 
((Xty)). The latter was not seen from Ventres' uniform thickness 
solution although the replacement of 6 in the kernel function 
was suggested, and the expansions, equations (33) and (37) are 
required to render the two integral equations (38) and (39) in 
conventional kernel function form. 

The expansion parameter in the above analysis has been chosen 
to be (1 - £ ). An alternate series expansion can be worked out 

using the parameter . The former expansion involves 

the detailed boundary layer thickness variation, while the latter 

only coricern^with the trailing edge boundary layer thickness 

variation. However, P in Equations (22) and (37) remains the 

1 

same for both choices of the expansion parameter. Although the 
latter expansion is probably formally preferable since a uniform 
series convergence is expected, the numerical results shown below 
have been carried out using the former. 
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IV. PHESENTATION OF COMPUTATIONAL RESULTS 

Tho shear flow lifting surface theory derived above has been used 
to calculate the pressure load on a two dimensional flat plate air- 
foil of finite chord in steady shear flow. Both Ventres* uniform 
thickness results and the present variable thickness results are 
presented for comparison. The familiar collocation method was used. 

The pressure was expanded as a linear combination of selected modal 
functions and the npKash was matched at all the collocation points 
The assumed modal functions, »^Cl-x/^)/Cx/^) for n = 0, 1, 2,...., 

are equivalent to but slightly different from those used by Ventres 
Tho details can be found in his paper and are not elaborated here. Shovm 
below are (a) the pressure distributions, (b) the total lift coefficients, 
and (c) the corresponding center of pressure travel for a turbulent 
shear flow model. In addition^ a lift coefficent resuic for a laminar 
shear flow model is given. For the turbulent shear flow case, the 
boundary layer was assumed to grow according to t$Cx)/6 = 

I B 3X C 

For the laminar shear flow case, it is assumed that 5(x)/6 =» 

'• max 

Ca) Pressure Distributions - Turbulent Boundary Layer 

Figures 4 and S show typical pressure distributions along the 
chord for two different ratios of the trailing edge boundary layer 
thickness to the chord length. The exponent N in the velocity 
power law is chosen as 7. It is seen that the variable thickness 
curve fits in between the potential flow (top curve) and the uniform 
shear (lower curve) results. Tho thicker the boundary layer is. 
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the larger the shear layer effect is on the pressure load. This 
is atainly due to the augmentation of the momentum thickness of 
the boundary layer which reduces the momentum transferred to the 
perturbed surface. Similar results for the N « 11 case have 
been calculated and show less shear effect because the N « 11 
case is closer to the potential flow CN » <^) than N « 7 case 
as far as the mean flow velocity distribution is concerned. 

Cb) Lift Coefficients « Turbulent Boundary Layer 

The normalized lift coefficients are plotted in Figures 6 

and 7 for N = 7 and N » 11 cases. Figure 6 indicates a 5\ 

increase in the total lift for the variable 6 compared to the 

uniform 6 case corresponding to • .1 and N » 7. Figure 

o 

7 shows less change as expected. The reductir>n of lift coeffici- 
ents is a natural result of the pressure decrease due to the 
existence of the shear layer. 

Cc) Center of Pressure Positions - Turbulent Boundary Layer 

Figures 8 and 9 demonstrate the shear layer effect on the 
center of pressure location. It is seen that the thicker boundary 
layers result in further backward movement of the center of 
pressure behind the quarter chord. This is caused by the more sig- 
nificant pressure reduction near the loading edge than that near 
the trailing edge for the shear flow. 

(d) Lift Coefficient for A Laminar Boundary Layer Flow 

Figure 10 shows the calculated lift coefficient for a 
laminar flow. The choice of N= 2.1 simulates the Blasius* 
laminar boundary layer velocity profile in an approximate way 
as shown in Figure 11. The choice of N = 2.1 instead of N® 2 is 
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because the kernel function is singular at N = 2. Due to this 
singularity the accuracy of the laminar result is somewhat suspect 
although Figi're lo shows the same tendency as the turbulent result. 
For N ■ 2 the present theory simply predicts zero lift. This is 
because the kernel function tends to infinity everywl.ere as N -*■ 2 
and hence one must have zero wall pressure to balance the finite 
upwash on the left hand side of the equation (38). This fact 
points out the limit of validity of the present shear flow model 
tor (nearly) laminar mean flow. Thus, one ought to include the 
viscosity effect, which is neglected in our shear flow model , in 
the analysis for the laminar boundary layer or> equivalently* postu- 
late a finite wall velocity according to Lighthill^^^ . 
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V. CCNCLUSION 

A nethod of calculating the pressure distribution cn lifting 
surfaces in incompressible shear flows with slowly varying boundary 
layer thickness has been developed < The assumed mean flow velocity 

I ' 

profile* u/y ■ (i/j(x,y)) is a good approximation for turbulent 
boundary layers at high Reynold's numbers* although it predicts an in- 
finite wall shear stress. An interesting point of using this velocity 
profile is the assumed zero mean flow velocity at the wall which avoids 
the problem of postulating a finite wall velocity as many authors have 
done. In any case the boundary layer reduces the pressure load com- 
pared to the potential flow* but the variable thickness shear flow 
model reduces it less for equal maximum shear layer thickness. 

The present theory for the variable thickness problem has been 
worked out in detail for the steady, incompressible case. Its extension 
to three dimensional* unsteady compressible flows seems to be workable. 
However* a question does arise as to whether the assumed pressure mode 
functions for supersonic flows should satisfy the kutta condition at 
the trailing edge. This is not clear because the flow is subsonic near 
the solid surface and is supersonic away from the surface. 
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AFPEI'DIX 


STEADY AIRFOIL MOTION IN A SHEAR LAYER OF VARIABLE THICKNESS 
We first consider for simplicity two-dimensional flow; it will 
be clear that the basic result hhlds in three -dimensions. Beginning 
with the equations of a fully viscous flow (Havier-Stokes equations) , 


>u ♦ 3w ■ 0 continuity (1) 

3x dz 


u^+w^3 - ) 3p 

3x 3z V Sx 


V / 3^u 


3^U 


streamwise 
X - component of 
momentum balance (2) 


u ^ + w ^ » - / J./ N \ transverse 

3x 32 ' j£ ' 3x^ Sz*^ / z - component of 

momentum balance (3) , 

We now construct an order of magnitude analysis of the various terms. 

For this purpose several scale factors are Introduced. The boundary 
layer thickness* S , is a characteristic length of the mean (shear) 
flow in the z direction; l is r characteristic leng; h in the x direct- 
ion associated with the growth of the boundary layer thickness; and a is 
a characteristic length in the x direction associated with the variation 
of airfoil upwash or angle of attack. First consider the mean flow with no 
iqtwash or angle of attack perturbation. We assume that 


3 1/ 

3x 


_ 1/ 
3z 




Ot) 
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From continuity » 


0 (u/j ^ +0 ( ) - 0 + w -w u 


(S) 


From X - momentum. 


( f ) ( r )” (^'0 


Using C^), the two terms on the left hand side of ((>) are seen to be 
of the same order. On the right hand side, the second term may be 
neglected compared to the last and the latter we require to balance with 
the left hand side 

Q / 0 / vu \ 


Thus 


6/j^ 0 ( Rey 


• 1/2 


(7) 


where the Reynolds number is given by 

Re s u4 

V 

From z - momentum. 


0 (») .0 ( 1 /^) | E *0 


CS) 


Using (5) and (7) , we conclude from C&) that 



REPR«-I)i;cr'''-T7 OR THE 
ORUm^AL X3 POOR 


- 0/p)H' ® (^) 

•nd 

- (ff) 


Thus 


3E« 

37 . 


3E 

3x 


and the aean pressure is essentially constant through the boundary layer. 
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These results are well'-knovm« of course, and due to Prandtl . We 
shall use a similar approach for the perturbation equations where 
there is a non-zero upwash or angle of attack. 

Let 

u (x,z) u(x»z) * fi(x,z) 

w (x,z) - w(x,z) + ft(x,z) (9) 

P Cx,z) = pCx»z> + ^Cx,z) 


where - denotes mean flow and “ perturbation. From our previous 
analysis ue note that FCx,zX ■ FCkX* Substituting C9) into Cl)> 
C2)^ C3) and subtracting out the mean flow equations one obtains the 
perturbation equations. These are 


sa ■» 3^ »» 0 

3x 3z 


continuity CJO) 
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v/32^i ♦ 32 q \ 


u 3ft ♦ 0 3w ♦ w 3ft » ft 3w 
3X 3X 3z 3z 


x>momentum (11) 


z -momentum (12) 

We as Slime that 

iC ) '^ ( )/, ; ; uQ-v C)/> 

32 3X “ 3X 

and 5 << t 

«/a 0 CD (13) 

a « t 



Other Interesting cases could be considered; for example, &/ » 1 
or << 1. However, we shall not pursue these here. 

From continuity, 

0 ^ft^ = 0-*-ft'u<i (14) 

Prom X - momentum. 
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In the above, we have used (5) and CI4) to replace w, ft and 
indicated by an arrow negligibly small terms. By the assumption of 
'v 0(1), the last two terms on the RHS of (IS) are of the same 
order. Let us compare the latter of these to the (remaining) terms 
on the LHS. Recall from (7) that 

Jt 

ut 

Thus 

vQ 'b ihl 

T 

which is negligible compared to LHS of (15). Hence, x-momentum equation 
siiq>lifies to 

u ^ ♦ ft - a/ ) ^ U6) 

Note that if a i , then one must include viscous terms and hence 
the basic validity of a shear flow, as opposed to a fully viscous 
model is dependent upon a « i. Formally tld? need not be true for 
turbulent flow where [22] 

UA 

However, then the question arises as to what viscosity coefficient 
one should use^®^. 


Prom momentum. 
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$ 


0 

V 


SI) * 



0 



- '*'.1 II * " 


C17) 


where we have used (S)* (7) and (14) to simplify terms in (17) and 
indicated by an arrow those which are negligibly small. Thus the 
z-nomentum equation becomes 


u 3ft - - (1/ ) ^ U8) 

3x ^ 3z 

Our final equations are then (10). (16) and (18), Using these we may 
obtain a single equation for p. 

V2p - 2 3U 3p - ^ 0 (19) 

3z 3z 

U ^ M 

where the last term may be neglected consistent with our previously 
announced assumptions on the various length scales. Note that we 
have dropped the and replaced u by U. 
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FIG. 6* LIFT COEFFIENT vs 8max/C FOR BOTH UNIFORM 
AND VARIABLE BOUNDARY LAYER THICKNESS 
OF THE CASE N * 7 




FK3. 7 . LIFT COEFFICIENT vs Smox/C FOR BOTH UNIFORM 
AND VARIABLE BOUNDARY LAYER THICKNESS 
OF THE CASE N= II 
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FIG. 10 -LIFT COEFFICIENT vs 8max/C FOR BOTH UNIFORM AND 
VARIABLE BOUNDARY LAYER THICKNESS OF THE 
CASE N = 2.1 







